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Abstract 

The Laplace transform is a powerful and versatile concept with broad applications in science and engineering. It 
allows a general means to represent and solve differencial equations, while being related to other important mathematical 
concepts such as power series and the Fourier transform. In this text, an introduction to the Laplace transform is 
provided, with emphasis on graphic representations, also including examples of differential equation solutions related 
to pendulum oscillations and linear electric circuits. 


‘II tempo rimane abbastanza a lungo per coloro che lo 
utilizzano.’ 

Leonardo da Vinci. 


1 Introduction 

Scientific modeling is, to a great extent, founded on mod¬ 
eling (e.g. pQ) the phenomenon of interest in terms of dif¬ 
ferential equations, typically followed by finding or simu¬ 
lation their solution. Examples of this approach abound 
in physics, acoustics, mechanics, electricity, magnetism, 
thermal sciences, among many other areas. Therefore, it 
would be great if one could find a way to a more general 
approach to handling and solving differential equations, 
which would otherwise require specific integration tech¬ 
niques. 

The good news is that such an approach exists: the 
Laplace transform (e.g. By transform¬ 

ing functions of a real variable t into a complex variable 
(e.g. 0) s considering a basis of complex exponentials 
e ~ st , this transform is capable of mapping equations and 
derivatives and integrals from linear differential equation 
into polynomials on s, therefore typically facilitating the 
respective solution. 

The one-sided Laplace transform of a real or complex¬ 
valued function /(£) of a real variable t > 0 is defined 
as: 


pOO 

C{fm*)=F{s)= f(t)e~ st dt. (1) 
Jo 



Figure 1: Some of the complex functions r](s) = e~ at e~ lujt defined 
by some respective values s = cr + iuj in the complex plane. The 
functions in blue are respective to cr = 0, corresponding to the same 
basis of the Fourier transform. The functions in green have cr < 0, 
therefore exhibiting exponentially decaying behavior. The functions 
represented in red correspond to cr > 0 and, consequently, increase 
exponentially. Observe the symmetric (‘mirrored’) organization be¬ 
tween the functions for uj > 0 and respectives —uj counterparts. The 
functions rj(s) are not shown to scale for better visualization’s sake. 

So, the Laplace transform maps the function /(£), t > 
0, into a complex function F(s) of a complex variable 
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s = cr + iw, where a and uj are real values. Keep in mind 
that this definition disregards the negative portion of the 
time domain of f(t) even if it exists, which is related to 
the causality intrinsic to the one-sided Laplace transform. 

Interesting insights about the essence of the Laplace 
transform can be inferred by considering the complex 
function r)(s) = e~ st of the complex variable s. Figure 
[l] illustrates this function for several values of s in the 
complex plane. 

Only the complex semi-plane defined by a > 0 is often 
considered in the Laplace transform, as several functions 
f{t) converge in this region. Each instance of cr and uj 
defines a respective complex function rj(s) = = 

e -at e iujt , represented graphically in Figure [I] For a > 0, 
we have that rj(s) decays steadily to zero for any value of 
cj, while oscillating according to uj. The limit situation 
where cr 0 is characterized by maintained oscillations. 

Because the term e~ at of rj(s) decays markedly with cr, 
the magnitude of the Laplace transform values obtained 
for increasing cr becomes smaller and smaller, in an expo¬ 
nential way, as illustrated in Figure [2j which shows rj(s) 
for cr = 0.1 and uj = 1. The decaying exponential envelope 
is shown in orange. 


This is valid only if Re(s) > 0 (i.e. a > 0). It follows 
that the Laplace transform of a constant function f(t) = c 
is c/s. We can also write 



1 

cr + iuj 


CT — iuj CT — iuj 

(cr + iuj)(a — iuj) ’ cr 2 + uj 2 


so that 


Re(F(s)) = 
Im(F(s)) = - 


cr 

cr 2 + uj 2 
UJ 

cr 2 + uj 2 


Figure [3] illustrates this Laplace transform graphically. 
The unit step function allows us to express the impor¬ 
tant fact that 


C{fms)=C{f(t)u(t)}(s). 

Another frequently used function is the causal expo¬ 
nential f(t ) = e~ at u(t), which has its Laplace transform 
calculated as 



t 


Figure 2: The exponentially decaying oscillating functions are the 
basis of the Laplace transform for cr > 0, where convergence is often 
observed. The orange curve corresponds to the exponential decaying 
function e~ at . 

A sufficient condition for the Laplace transform of f(t) 
to exist is that this function has exponential order , mean¬ 
ing that 

| e -VW| < ce s (s-«y 

for a real constant c and a real value a > 0 so that 
Re(s) = cr > a. 

As an example, the Laplace transform of the unit step 
function ( u(t ) = 1, t > 0 and 0 otherwise) has as Laplace 
transform: 



L(s) = / e- at e~ st dt = 

Js =0 

for Re(s) < —a. 

Table [l] lists several Laplace transforms of often consid¬ 
ered functions, including the respective regions of conver¬ 
gence. 

As it can be observed, the Laplace transforms tends 
to exhibit a ratio form, with the complex variable s in 
the denominator (i.e. l/(s + a)). This is also reflected in 
that visualizations of Laplace transforms tend to have an 
accentuated decay for increasing values of cr (i.e. along the 
real axis, as in Figure |3|, though with varying number of 
peaks (e.g. the cosine and sine functions have two peaks 
in the respective Laplace transforms). Perhaps for this 
reason, Laplace transforms are not frequently visualized, 
though the position of roots and poles is often presented 
graphically in the complex plane. 

The two-sided Laplace transform is defined as 

/ oo 

f{t)e~ st dt. 

-oo 

Observe that 

C{f{t)u(t)}{s) = C{f(t)}(s). 

The Fourier transform , a non-causal transform that is 
related to the steady-state regime, is related to the two- 
sided Laplace transform in the sense that J 7 {f(t)} = 
F(s) |cr=o. We also have that 


(s + a) 


0 (s—a)t 
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/ oo 

f(t)e- at e- M dt = T{e- at m}. 

-OO 

The Laplace transform is also related to other trans¬ 
forms (e.g. the Mellin and ^-transform) and mathemat¬ 
ical concepts (e-g. 0 ) such as statistical moments and 
the power series of a discrete function h[n\, namely 
h[ n ] xn ( n is a positive integer and x is a real vari¬ 
able), corresponding to the limit situation when n —> t. 


In an similar way, considering that sin(u)t) = (e luJ — 
e~ luJ )/{2i ), we have that 

_ r • / m w 

We can also apply the complex shifting property to eas¬ 
ily obtain the Laplace transform of e~ at cos(ujt ) as 


C {e at cos(co’t)} 


w 

(s + a ) 2 + uo 2 


Table[2]lists several properties of the Laplace transform. 


2 Properties of the Laplace Trans¬ 
form 

The Laplace transform has several properties that are par¬ 
ticularly useful. For instance, the complex shifting prop¬ 
erty says that 

C [e at m] = F(s - a), (3) 


3 The Inverse Laplace Transform 

The inverse of the Laplace transform, namely £ _1 {F(s)}, 
so that £ _1 {£{/(£)}} = /(£), can be defined as 

1 />7+iT 

m = ins)} (t) = — lim / F(8)e«d8, 

Am oo J 7 _ iT 


where a is any complex number. 

We can combine Equations [ 2 ] and [3] to derive 


C{e at } 


1 


s — a 


In the particular case that a = icj, we have 


C{e iuJt } 


1 


s — iuj 


Now, we have that the cosine function can be expressed 
as cos(u;t ) = ( e luJt + e~ luJt )/2 , so that 


where 7 is a fixed real part of s, so that the integration 
proceeds along the respective defined straight line (paral¬ 
lel to the imaginary axis). 

Because this integral can be difficult to be calculated, 
the inverse Laplace transform is often obtained by orga¬ 
nizing F(s) in a convenient way, so that tables of the 
Laplace transforms (such as Table [T]) can be applied in 
order to to derive the respective inverse. Examples of 
this procedure will be presented in the remainder of this 
text. 


C {cos(oot)} =-j- 

S — VjJ 

1 s T iuj T s — iuj 

2 (s + iuj)(s — iuf) 


■ VjJ 

s 


S 2 T id 2 * 
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4 Ordinary Differential Equations 

Let’s consider a generic second order ordinary differential 
equation (e.g. 0) 























a 


im 

dt 2 


+ b 


df(t) 

dt 


+ c = g(t), 


where the function g(t) represents an external effect 
driving the dynamics of the modeled system. When 
g(t) = 0, we say that the equation is homogeneous. 

The Laplace transform of this equation yields 


The Laplace transform-based solution of this equation 
follows immediately as 


sF(s) - /(0) - aF(s) = 0 => F(s) = -tM. =*. 

s — a 

=> fit) = /(0)e“*. (5) 


s 2 F(s) - sf( 0 ) - /( 0 ) 


+ b[sF(s)-f(0)}+c = G(s). 


Thus, we can write 


6 The Harmonic Pendulum 

The characterization of small oscillations of a pendulum 
can be approximated by the linearized (small oscillations) 
differential equation (e.g. 0). 


_ (as + b)f(0 ) + af(0 ) ^ G(s) 

as 2 + bs + c as 2 + bs + c 

and it follows that 


F(s) 


P(s) 

Q(s) 


A ( s ) , G(s) 
Q(s) Q(s)' 


(4) 


Q(s) is called the characteristic function associated to 
the original differential equation. The term can be 
understood as the Laplace transform of the homogeneous 
equation (i.e. excluding the driving function). 

Given a ratio R(s) = of two complex function, the 
roots of Q(s) are called the respective poles of i?, while 
the roots of M(s) are called simply the roots of R. 

Let Q(s) in Equation [4] have pi all mutually distinct 
poles pi. In this case, it can be shown that 


P(s) _ _ P(s) _ = 

Q(s) (s - pi ) (s -p 2 )...(s- p M ) 

— Cl ° 2 + ° M 
s-pi s-p 2 s — p M ’ 

where the constants q can be obtained (e.g. 0) as 




{s- Pi )P(s) 
s^Pi Q(s ) 


This procedure, called partial fraction expansion , will 
be illustrated in several of the following sections. Note 
that a different procedure applies when we have two or 
more identical poles. 


4>{t) + c<t>{t) = g(t), 

where c = g/L, and g is the gravity acceleration and L 
the length of the pendulum rod. 

We can write 


s 2 $(s) — s4>{ 0) — 0(0) 


+ c<f>(s) =G(s) => 


s 2< f>(s) + c$(s) = s0( 0) + 0(0) + G(s) => 

_ S( t ) i^) + ^( 0 ) | G( s ) 

[ ’ s 2 + c + s 2 + c 

_ Pjs) _ Ajs) G(s) 

Q(s) Q(s ) Q(s)' 


Let’s consider that g(t) = 0, so that 


<r>M _ ^( s ) = 5 0 (O) + 0 ( 0 ) 

U Q(s) s 2 + c 


Assuming 0(0) = 0 and 0(0) = 0o, we get 


$(S) s 2 + wg ’ 
where cjo = y/c = y/g/L. 

The solution of the original differential equation can 
now be found by taking the inverse Laplace transform 
(through Table [l]), yielding 


(j)(t) = 00 COS(cuot) = 00 COS 



7 The Driven Pendulum 


5 Exponential Growth 

One of the simplest ordinary differential equations is de¬ 
fined by the hypothesis that the rate of growth of a 
given function f(t) is proportional to its current value 
(e.g. Malthus’s population model), i.e. 

j , m = am 


Another interesting situation involving a linear approxi¬ 
mation of a pendulum (small oscillations) is when an ex¬ 
ternal force g(t) is applied to drive the movement. This 
can be expressed as 

0(f) +c0(t) =g(f). 

As before, let’s consider that 0(0) = 0 and 0(0) = 0o, 
yielding 
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\ = , C(s) _ ^( s ) , G( s ) 

s 2 + c s 2 +c Q(s) Q{s)' 

We already obtained the inverse of A(s)/Q(s ), given as 
in Equation |6j which we will call <pA (t) = 0o cos(ujot). For 
g(t) = asin(ujt ), we have 

G(s) = THZ+ = cluj _ = 

Q(s) s 2 + c (s 2 + uj 2 )(s 2 + c) 
auj 

(s - Pi)(s - p 2 )(s - p 3 )(s - Pi) 

— Cl _|_ ° 2 q_ c 3 C 4 

S - Pi S -p 2 S -p 3 S -p 4 

Partial fraction expansion yiedls 

G(s) CLUJ 

Q(s ) (s + iuj)(s — iuj)(s + iy/c)(s — iy/c) 

Since (s 2 -\-uj 2 ) = (s + zu;)(s — icj), we immediately find 
that pi = — iuj and P 2 = ++j. The other two poles can 
be found by solving the quadratic equation s 2 + c = 0 as 
P 3 = —iy/c an d Pa = i yfc- Therefore, it follows that 


<^(t) = Cl e~ iuJt + c 2 e iuJt + c 3 e~ iV ~ ct + c 4 e‘ i ^~ ct . 

Recall that ujo = y/c. The constants c* can be calcu¬ 
lated as 


Thus, we find that the driven pendulum will move as a 
superimposition of two oscillations with respective angu¬ 
lar frequencies ujq and uj. Figure [4] illustrates the resulting 
movement 0(£) assuming 0o = 7r/200, g = 10, L = 0.1 (so 
that cjq = 10 rad/s), a = 0.2, and uj = Ylradjs. 



t 

Figure 4: The oscillations of the small movement pendulum driven 
by a sinusoidal function, assuming </>o = 7r/200, g = 10, L = 0.1 (so 
that ujq = 10 rad/s), a = 0.2, and lj = Ylradjs. 


8 The Damped Pendulum 

Let’s now consider that the pendulum movement is sub¬ 
jected to some damping (e.g. air resistance), which can be 
modeled as 




Let’s make a = f and (3 = and write 


ci = 


2 uJq — uj 2 


c 2 = -- 


2 UJq — UJ 2 


-Cl 


a 

°3 = x 


IUJ 


2 uj 0 {uj? 


a 

C A = -X 


IUJ 


2 cj 0 (^o ~ ^ 2 ) 


o-^ 2 ) 
= -c 3 . 


Now, we can write the solution respective to the term 

G ( s ) as 
Q(S) 


cjjsit) = cie-^ - ae iujt + c 3 e" <Wot - c 3 e^ ot = 


= Ci [e~ iujt - e iujt ] + c 3 [e 


—i ujot 


— e 


iujotl _ 


% a 


2 K-u 2 ) L 


( e -iu,t _ e iu,t\ + f e -iuot _ e iuot\ 

W 0 

a cluj 

—o- ^sin(ujt) H-7—9-y -sin(uJot ). 

UJq — UJ 2 UJo{UJq — UJ 2 ) 


So, the full solution can be expressed as 


0(£) + acj){t) + 00(£) = 0 =+ 
s 2< f>(s) — s0(O) — 0(0) + a$(s) + f3 [s$(s) — 0(0)] = 0 =+ 
s 2< f>(,s) + f3s<fr(s) + a$(s) = 50(0) + 0(0) + /30( 0). 

Let’s consider that 0(0) = 0 and 0(0) = 0o, so that 


*(«) = 


(s + 0)0Q 
s 2 + /3s + a 


If 0 2 > 4a, we have two distinct real poles pi = 
(-0 + ^/0 2 - 4o) /2 and p 2 = (^-0 - f3 2 - 4a^ /2, 

and it follows that 


®(s) = 0o- 


s + /3 


(s -Pi)(s -p 2 ) 
and we find 


— 00 


ci 


C2 


+ 

s-Pi S-P 2 


4>(t) = 4> 0 [cie Plt + c 2 e P2t ] , 
with pi < 0, P 2 < 0, and 


= 00 cos(uJot) + 


CLUJ 


cj 0 (^o ~ ^ 2 ) 


0 (t) = 0 T (t) + 0 5 (t) = 

sin(uJot) H-— 7zsin(ujt). 


UJn — 


Cl 


lim 


S + 0 
5 -P 2 


Pi+j9 
Pi P 2 


10+ y/ f3 2 — 4a 
2 ^/0 2 — 4a 
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r s + /3 P 2 + /3 1 /3 - y/fi 2 - 4a 

s^p 2 S~p 1 P2-P1 2 y^ 2 _ 4a 

This parameter configuration is often called over 
damped , with the oscillations decreasing to zero exponen¬ 
tially, without oscillations, as implied by Equation [8] 

Two other situations could be considered, namely: (i) 

/ 3 2 = 4a (two identical real roots) and /3 2 < 4a (two com¬ 
plex roots), which correspond respectively to the critically 
damped and under damped configurations. In the former 
case, the pendulum returns to the equilibrium as fast as 
possible, without oscillating and never reaching the equi¬ 
librium position. In the latter situation, the pendulum 
oscillates with decreasing amplitude as it returns to the 
equilibrium. Finding the respective solutions is left as an 
exercise. 


= sV(s) -v c (0) = 

L(s + a + iu>)(s + a — iuS) 

_ sV(s) - vc(0) 

L [(s + a) 2 + uj 2 ] ' 

For v(t) = Au(t), we have that V(s) = A/s, and 


I(s) 


A - t>c(0) 

L [(s + a) 2 + uj 2 } ’ 


Thus we find the solution 


i(t) = -— 

(jj-L 

Figure [5] illustrates this solution for R = O.Icj, L = 
1 mH, C = lOraF, fc(0) = 0.5V, so that i(0) = 0.5 /R = 
5A. 


9 RLC Series Circuit 

Laplace transforms are widely used to solve linear circuits 
(e.g. [10]) involving resistors (R), capacitors (C) and in¬ 
ductors (L). In this section we apply the Laplace trans¬ 
form to solve a circuit composed of a voltage source v(t) 
in series with a resistor, a capacitor and an inductor, as 
well as a switch that can be in open or closed position. 
The circuit equation in the time domain is 

If 1 d 

Ri(t) + — J i(t)dt + v c {0) + L — i(t) = v(t). 



t 

Figure 5: The solution i(t) of the RLC series example for R = 0.1a;, 
L = ImH, C = 10 mF, vc( 0) = 0.5V, so that i(0) : 0.5 /R = 
5 A. Observe that the current i(t) undergoes decreasing oscillations, 
converging steadily to zero. 


which yields 


RI{s) + ^— + + l [s I(s) - i(0)] = V(s). 

C s s 

Assuming that i(0) =0 (.e. the switch is open until t = 
0, allowing the capacitor to maintain its initial condition 
vc(t) before that moment), we have 


R 4~ 


i 

Jc 


T sL 


I(s) = V(s) - 


Vc( 0) 

s 


10 Transient and Steady-State 

Let’s consider the previous RLC circuit, but now con¬ 
nected to a sinusoidal voltage source with fixed frequency, 
as it is often the case in electrical analysis. More specif¬ 
ically, onsider that g(t) = v(t) = sin{ujt ), implying 

V(*) = TO- 

We have from the previous section that 


Hence 


1 ( 8 ) = 


sV(s) - V C (0 ) 


L (s 


T s 


Tc) 


The poles of the characteristic equation can be obtained 


as 


s = 




—a dz iuj. 


Assuming > (^;) 2 , it follows that 


I(s) 


^c(q) 

L [(s + a) 2 + uj 2 ] 


which can be developed as 


t(\- SU} ~ s2v c(0) - uj 2 v c (0) _ 
“ L(s 2 +u 2 )[(s + a) 2 +uj 2 ] ~ 
1 suj — s 2 vc{ 0) — uj 2 vc(0) 

L (s + iuj)(s — iuS) [(s + a) + iu\ [(5 + a) — icj] 

Applying partial fraction expansion, we obtain 
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SUJ — 8 2 t’c(0) — U 2 Vc(0) 

(s + iw)(s — iuS) [(s + a) + iuo\ [(5 + a) — iuj\ 

_ _Ci_ _C2_ c 3 | c 4 

s + iuj s — iuj (s + a) + iuj (s + a) — iuo 
The respective solution is given as 


*(i) = y \c\e~ iut + c 2 e iujt + c z e~ at e~ iut + c A e~ at e iujt ] , 

±J 

For a sufficiently large time t —> oc, the transient stage 
will have mostly vanished and the circuit will be in its 
steady-state regime, implying that 

i s teady(t) = ^[c 1 e- Mt +c 2 e Mt }. 

So, in the steady state, the current is a sinusoidal func¬ 
tion with the same frequency as the voltage source. In this 
situation, the effect of the original charge in the capaci¬ 
tor will have almost completely vanished. This behavior 
allows us to consider the solution of steady-state circuits 
by using phasors (e.g. ED- 


11 Concluding Remarks 

In this text, we provided an introduction to the Laplace 
transform and applications to solving differential equa¬ 
tions. It should be kept in mind that this transform can 
have many other interesting applications, such as solving 
some types of partial differential equations, circuit net¬ 
works, as well as being applied to many other problems. 

In particular, the so-called lumped approach (e.g. [12]) 
allows the circuit analysis methodology illustrated here to 
be applied to other areas such as as mechanics, acoustics, 
thermal sciences, magnetism, among many other possibil¬ 
ities. 
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CDTs intend to be a halfway point between 
a formal scientific article and a dissemination 
text in the sense that they: (i) explain and 
illustrate concepts in a more informal, graphical 
and accessible way than the typical scientific 
article; and, at the same time, (ii) provide more 
in-depth mathematical developments than a more 
traditional dissemination work. 

It is hoped that CDTs can also provide integration 
and new insights and analogies concerning the 
reported concepts and methods. We hope these 
characteristics will contribute to making CDTs 
interesting both to beginners as well as to more 
senior researchers. 

Though CDTs are intended primarily for those 
who have some preliminary experience in the 
covered concepts, they can also be useful as 
summary of main topics and concepts to be learnt 
by other readers interested in the respective CDT 
theme. 

Each CDT focuses on a few interrelated concepts. 
Though attempting to be relatively self-contained, 
CDTs also aim at being relatively short. Links 
to related material are provided in order to 
complement the covered subjects. 

The complete set of CDTs can be found 
at: https : //www.researchgate.net/pro j ect/ 
Costas-Didactic-Texts-CDTs, 
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